o 
o 
o 

(N 



(N 



> 



Nonequilibrium perturbation theory for spin-^ fields 

I. D. Lawrie and D. B. McKernan* 
Department of Physics and Astronomy, The University of Leeds, Leeds LS2 9JT, England 

(July 12, 2000) 

A partial resummation of perturbation theory is described for field theories containing spin-i 
particles in states that may be far from thermal equilibrium. This allows the nonequilibrium state 
to be characterized in terms of quasiparticles that approximate its true elementary excitations. In 
particular, the quasiparticles have dispersion relations that differ from those of free particles, finite 
thermal widths and occupation numbers which, in contrast to those of standard perturbation theory 
evolve with the changing nonequilibrium environment. A description of this kind is essential for es- 
timating the evolution of the system over extended periods of time. In contrast to the corresponding 
description of scalar particles, the structure of nonequilibrium fermion propagators exhibits features 
which have no counterpart in the equilibrium theory. 
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I. INTRODUCTION 



Many interesting physical problems, arising, for example, in the study of the early universe and in relativistic 
heavy-ion collisions |^,^ require an understanding of the evolution with time of highly-excited states of a quantum 
field theory. The properties of high-temperature states in thermal equilibrium have been studied for a long time and 
much is known about them (see, for example ||^). For the most part, attempts to study the nonequilibrium properties 
of systems that evolve with time have been based on the assumption that this evolution can be adequately represented 
I . as a sequence of near-equilibrium states. While such an assumption may be justified in some cases, our understanding 
' of the true nonequilibrium dynamics is at present very incomplete. One route towards a more complete understanding 
, that has been pursued by several groups is through the study of the N ~^ oo limit of A^-component scalar field theories 
' (see and references cited in these papers). The major advantage of this limit is that it is a gaussian field theory 
for which the path integral can be evaluated exactly, yielding closed-form evolution equations suitable for numerical 
■ solution. However, it is also a theory devoid of scattering processes, which in general lead to important dissipative 
Qh! effects (and the same is true of the related Hartree and one-loop approximations). Moreover, it seems to be extremely 
^ difficult to extend such calculations beyond leading order in 1/N, which is essential for approaching any description 
^ ' of more realistic systems. In this paper, we focus on the most obvious alternative of extracting as much information 
, as possible from perturbation theory. 

A serious limitation of standard perturbation theory is that, being an expansion about a non-interacting theory, it 
too is devoid of scattering at lowest order. This means, in particular, that the occupation numbers of single-particle 
' modes which appear in propagators are fixed at some initial values and do not reflect the evolving nonequilibrium 
H ' state. Low-order calculations therefore become essentially meaningless unless they are restricted to time intervals 
. . . 1 much shorter than a typical relaxation time. To improve this situation, one should reformulate perturbation theory 
so as to describe the nonequilibrium state in terms of its own quasiparticle excitations. These excitations have a 
nonzero thermal width, which in part also characterizes the rate of relaxation of their occupation numbers in response 
to a changing environment. To put this idea into practice, it is necessary to construct a lowest-order approximation 
to the interacting theory in which at least some of the dissipative effects of interactions are resummed. Methods for 
achieving this in the case of both real and complex scalar theories have been described in |p^-p^ and incorporated in 
a comprehensive perturbative approach to the nonequilibrium dynamics of phase transitions in [ p3[ . The purpose of 
this paper is to investigate how the same idea might be implemented for spin-i fields. The nonequilibrium dynamics 
of spinor fields turns out to be quite complicated. In contrast to scalar fields, their propagators appear to have a 
structure that is not simply a time-dependent generalization of the one that applies in thermal equilibrium; it is 
sufficiently complicated that we have not been able to explore it in full generality. 

We begin in Section II by reviewing briefly the resummation of 2-point functions for real scalar fields. In Section III, 
we derive some general properties of the full spinor 2-point functions which serve as a guide to the construction of an 
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effective quasiparticle action, under the simplifying assumption that the latter will be CP-invariant. The quasiparticle 
action is constructed in Section IV in terms of several undetermined functions of time and spatial momentum that 
characterize the quasiparticle dispersion relation, thermal width and occupation numbers. These functions appear 
in a counterterm which is added to the free part of the action and subtracted from the interaction part, and will 
subsequently be determined self-consistently by requiring the counterterm to cancel part of the higher-order corrections 
to the self energy. The real- and imaginary-time quasiparticle propagators corresponding to this effective action are 
derived in Sections V and VI respectively and a self-consistent criterion for determining the quasiparticle parameters 
is implemented in the context of a simple model in Section VII. For illustrative purposes, we introduce supplementary 
approximations that allow them to be evaluated in closed form. These approximations correspond to a weakly 
interacting system close to equilibrium, and for this special situation we find, reassuringly, that the evolution of 
occupation numbers is described by a Boltzmann equation. Finally, in Section VIII, we summarize our principal 
conclusions and comment on their relation to some other approaches to non-equilibrium field theory. 



II. DISSIPATIVE PERTURBATION THEORY FOR SCALAR FIELDS 



Consider the usual Xcj)^ theory, defined by the Lagrangian density 



(2.1) 



and suppose that an initial state of thermal equilibrium with inverse temperature (3 is set up at time t — Q. In this 
model, the time-dependent mass m{t), which arises, for example, in the case of a scalar field theory in a Robertson- 
Walker spacetime, drives the subsequent state away from equilibrium. Then the closed-time-path formalism (developed 
for general time-dependent situations in EJ]) yields a path integral weighted by the action 



3) = 



(2.2) 



where the path- integration variables (jji, (j32 and 03 live on a closed contour C in the complex time plane. This contour 
runs along the real axis from i = to a final time tj, returns along the real axis to i = 0, and finally descends along 
the imaginary axis to t = —ip. Here, we have taken the limit — > 00. The Euclidean action £e (which uses m(0)) 
represents the initial density matrix. In this theory, there is a 3 x 3 matrix of 2-point functions Gab{x,x'), with 
a,b = 1, 2, 3, but our attention will focus mainly on the real-time functions, with a,b = 1,2. For the real-time part of 
the action, we wish to construct a lowest-order version /co(0ii ^2) = ~5 / di^x cjia'Dab'Pb, where, after a spatial Fourier 
transform, the differential operator T) is 



k'^ + m'^{t) 

-[df + P + m^{t)]. 



Mkit,dt 



(2.3) 



The counterterm i J dl^x (j)aMab4'b is added to /co and subtracted from the interaction /cint = ~ ^CO: so as to 
leave the whole theory unchanged. A choice of is a choice of the approximate theor y a bout which we perturb and 
is, of course, equivalent to a choice oiV. Subject to several constraints (discussed in [Q, and generalized below for 
spinors), the most general choice for T) is 



Pk{t) - iakit)] hk{t)dt + ^jk{t) + iakit)] ' 

-Ik {t)dt - i-Tfe (t) + iak {t)] [~dl - (3k (t) - iak{t)] . 



(2.4) 



where ak( t)^ (3 k(t) and 7fe(t) are real functions yet to be determined. Of course, the counterterm M can be read from 
(U) and (g]|). 

The 2x2 matrix of quasiparticle propagators gk{t,t') is the solution (subject to suitable boundary conditions) of 

Vkit,dt)gk{t,t') ^ gk{t,t')Vk{t',^tt') = -i5{t-t') . (2.5) 
Suppressing the spatial momentum k, this solution can be written in terms of a single complex function h{t,t') as 

9ab{t, t') = hb{t, t')e{t - t') + ha{t\ t)e{t' - t) , (2.6) 
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where hi = h and h2 = h*. The function h is 



h{t,t') = iexp 



Jit")dt"j {[l + Nit')]f(+\t)f -\t') + [-l + N*it')]f(-\t)f+\t')} , (2.7) 
We see that one of the undetermined functions, jk{t), can be 



with = [2n{t)]-^/^e^p(^T^Jo^{t")dt" 

interpreted as a quasiparticle width. The quasiparticle energy fifc(t) is a solution of 

liifc(t) 3f2|(t) „o,.. 1 



2r!fc(t) 417|(t) 



(2.8) 



FinaUy, the function Nk{t), which we hope to interpret in terms of time-dependent occupation numbers, is a solution 
of 



nk{t) 



[dt + -fk{t)] Nk{t) ^ 2tak{t) . 



(2.9) 



To give subs tanc e to this scheme, a prescription is needed for determining the three functions ak (t) , (3k (t) and jk {t) 
introduced in (2^). To this end, define the 2x2 self energy matrix T,k{t,t') by 



Gk{t,t')=gk{t,t')+i / dt"dt"'5fc(t,i")Sfc(t",t"')Gfc(i"',i')- 



(2.10) 



This self energy has contributions from the counterterm A4 and from loop diagrams: 

j:k{t,t')^Mk{t,dt)S{t-t') + ^'°°^{t,t'). (2.11) 

The general strategy is to optimize gk{t, t') as an approximation to the full two-point functions Gk{t, t') by arranging 
for A4 to cancel some part of I]'°°p. Clearly, since I]1°°p is non-local in time, only a partial cancellation can be 
achieved. Various prescriptions might be possible; perhaps the most obvious is the following. Express T,k{t,t') in 
terms of the average time t = ^{t + t') and the difference {t — t') and Fourier transform on [t — t'). The components 
of M.k{t,dt) contain at most one time derivative, so the self energy can be decomposed into contributions that are 
even and odd in the frequency: 



,(2) 



Generalized gap equations to be solved for ak{t), Pk{t) and 7fc(i) can now be obtained by requiring 



Afi^)(t-)^-4^) {mm 



(2.12) 

(2.13) 
(2.14) 



which amounts to an on-shell renormalization prescription. 

These gap equations provide exact implicit definitions of Q;fe(t), /3fc(i) and jk{t), but they cannot, of course, be 
exactly solved. If the perturbative expansions for E^^^ and s[,^^ are truncated at some finite order, one 
obtains concret e ex pressions for them in terms of the propagators gkit^t'). These truncated gap equations, together 
with equation (2.8) for the quasiparticle energy and (2.9) for the function Nk{t) form a closed system that one might 



try to solve numerically. It is to some extent illuminating to establish a connection with kinetic theory through 
some further approximations. Suppose that the gap equations are truncated at two-loop order - the lowest order 
that yields a nonzero quasiparticle width ^k{t)- Then, assuming sufficiently weak coupling and sufficiently slow time 
evolution, propagators inside the loop diagrams can be approximated by taking J*, ilk{t")dt" « Qk{t){t — t') and the 
limit 7fc(i) — > 0. Then, with quasiparticle occupation numbers nk{t) defined by 



Nk{t) 



l + 2nk{t) 



l-ijk{t)/2nk{t) 



(2.15) 



a time-derivative expansion of (2.9), yields the Boltzmann-like equation 



dtnk{t) 



A2 



32(27r)E 



d^fcid^/fc^d^fca^^ 



- - Qk) 5 (ki 



k) 



fllQ.2^3^k 

X [^1^2(1 + n3)(l + Uk) - (1 + ni){l + n2)n3nk] 



(2.16) 



In the following sections, we investigate how this resummation scheme might be extended to spin-i fields. 
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III. EXACT PROPERTIES OF SPINOR 2-POINT FUNCTIONS 



To be concrete, we consider a system defined by the Lagrangian density 

C ^Ip [j^'df, ~ m{t)]ijj + AC 



(3.1) 



where AC represents the coupling of the spinor ijj to other fields. For the purposes of this work, we suppose once more 
that the system is driven away from thermal equilibrium by the time-dependent mass m(t) (and possibly by other 
time-dependent parameters in AC), but remains spatially homogeneous. As in the scalar case, spinor field theory in 
a flat Robertson- Walker universe can be represented as a Minkowski-space theory with time-dependent mass. If the 
state at an initial time that we shall call t = is one of thermal equilibrium, then standard methods (described, for 
example, in Ref. serve to derive the generating functional 



Z{C,0= / l[[di>M exp[z/e(V,V')+*Se(^,^,^,C)] 

a=l 



(3.2) 



which generalizes that described in Section II for a scalar field. In this case, the source term is 

5c(?,V',?,0 = /*' dt \^At)Mt) +Mmi{t) +Ut)Mt) + i^2m2{t)] + f dT \;U^)Mr) (3.3) 



and we do not indicate explicitly the other fields that may appear in AC. If the initial state is characterized by a 
temperature (3~^ and chemical potential ^, then the path-integration variables at the ends of the closed time path 
obey the boundary conditions "01(0) — —e^^ipsiP) and V'i(O) — —e~^^ip^{f3), which are inherited by the Green's 
functions. 

As before, we are particularly concerned with the real-time 2-point functions 



■InZ 



for a, 6 = 1, 2. In terms of field operators, they are 

s^:l;\x,t;^,t')=[ 

\ {ij^{x,t)i,f,{x',t'))^ {T[^^{x,t)i^p{x',t')])^,. 



(3.4) 



(3.5) 



where a and (3 are spinor indices, while T and T denote time ordering and anti-time ordering respectively. In the 
presence of a chemical potential, the expectation values are given by 



Tr 



Tr 



(3.6) 



where H is the Hamiltonian at the initial time and N is the particle number. 

We hope to construct a perturbation theory in which the lowest-order propagators are partially resummed versions 
of these full 2-point functions, and begin by establishing some properties of the full functions that our approximate 
ones ought to share. Expecting that correlations should decay, very roughly as e^'^'*^* over large time intervals, we 
write the Wightman function S^p{x, t; x' , t'\ ^) = {ipa{x, t)ipp{x' , t'))i^i as 



S>p{x, t- x', t'; /i) = Ha^ix, t; x', t')0{t - t') + Ko.p{x, t; x\ t')e(t' - t) . 



Using ipfj — "0777/5; it is easy to see that 



and hence that 



S^^'ix, t; x', t'; fi) = -f"S>{x', t'; x, t; ^)7° 



IC{x, t; x' , t'; fi) = H{x' , t'; x, t; fi) , 



(3.7) 

(3.8) 
(3.9) 
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where, for any Dirac matrix, we define M = j^M^j^. It would be helpful if the second Wightman function 
S^p{x,t;x' ,t'; ^) = —{ipp{x.\t')'ijja{^,t))fj_ could be expressed in terms of the same matrix H{x,t;x',t';fj,). This 
can, in fact, be done in a CP-invariant theory. If the CP transformation of an operator A is implemented by a unitary 
operator Ucp, so that A'^^ = U^^pAUcp, then a CP-invariant theory has H'^^ = H and N*^^ = —N, and we see from 
that 

iCP 



{A),^{A^^)^,. (3.10) 



For a Dirac spinor, we have ip'~^^{x, t) — joCtp {—x, t), where C is the charge conjugation matrix and ^ indicates the 
transpose. It follows from this that 

S<p{x,t;x',t';fi) = [c-'joS>pi~x',t':-x,t;-fi)^°cY 

^ n{-x, t; -x',t'; fi)e{t - t') + ni-x', t'; -X, t; n)e{t' - t) , (3.11) 

where, for a matrix-valued function of the chemical potential, we define M{ii) — [C~^7°Af (— /x)7°C] It is simple 

to check that M{ii) — M{fi) = M{ii) and that M{ii) = M(/i). With these definitions, the matrix of real-time 2-point 
functions for a fermion in a CP-invariant theory can be expressed, after a spatial Fourier transformation, as 

5(-)(M';fc)=f'^^''''^^ ?(^'^'^)V(t-0+fS!*?'S (3-12) 

^ ' \n{t,t';k) nit,t';k)J ^ ' \nit',t;k) ■Hit',t;k)J ^ ^' ^ ' 

with 

n{t,t';k)= J d^x e-'''''n{x,t;0,t') . (3.13) 

Equivalently, defining n'''^\t,t' \k) = n{t,t';k) and H^'^'>{t,t';k) = n{t,t';k), we can write 

^^''^^ {t, t') = n^^\t, t'-k)e{t - t') + H^''\t', t; k)0{t' - t) . (3.14) 

We shall demand of our perturbative propagators that they have the structure shown here (as, indeed, do the 
propagators of standard perturbation theory). This does not mean that our resummation can be applied only in 
the context of a CP-invariant theory; it does mea n, th ough, that any CP-violating effects will not be resummed. It 
is worth pointing out that a relation similar to (3.11) can be obtained by assuming C invariance rather than CP 
invariance. This would be equ ally u sabl e, bu t phenomenologically a little more restrictive. 

The structure expressed by ( ^.12| ) or ( [3.14 ) implies two symmetries that will be useful to us. They are 

S^"-''\t,t'-k)^S^'"'\t',t-k) (3.15) 

and 

(S^'^\t,t';k) 5^''\t,i';fe)\ /5(22)(t',i;fc) S'^^^\t' ,t; k)\ 
\S^^^\t,t';k) S^^^\t,t':k)) V5(2i)(t',f;fc) S^^'\t' ,t; k) ) ' 



(3.16) 



The first of these generalizes to the full 3x3 matrix of real- and imaginary-time 2-point functions, but the second 
makes sense only for the real-time functions. 

Finally, we shall need two pieces of information concerning the values of these functions at equal times. The 
functions 5^^^-' and iS*^^^-* have unique values at t = t' , which implies 

n(t,t;k) ^nit,t;k) . (3.17) 

On the other hand, the time-ordered function S^-^^^ has a discontinuity a.t t = t' , which reproduces the equal-time 
anticommutator 



and this implies 



^^hm ^5^^^^^ {x, t; x', t') - ^,nm ^^^"^ t; x' , t') = ({^„((r, t), V^^(a;', t)}) ^ = j"^p5{x - x') , (3.18) 



n{t,t-k)-n{t,t-k) ^-1° . (3.19) 
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IV. CONSTRUCTION OF THE QUASIPARTICLE ACTION 



We wish to construct a lowest-order action 



(4.1) 



(2) 

that will serve as a starting point for our partially resummed perturbation theory. As before, 1^ is the quadratic part 
of the original closed-time-path action, while the term involving Mab is a counterterm which will be subtracted from 
the interaction part, so as to leave the whole theory unchanged. Specifying the form of ^Aab is equivalent to specifying 
the resulting differential operator "Dab- To begin, we construct the real-time components of 2?ab, with a, 6 = 1, 2. The 
unperturbed propagator matrix S'^"''-' (t, t') is a solution of the equations 



(4.2) 



(For economy of notation, we shall usually not indicate explicitly the dependence of these quantities on k and /i.) 
The form that might usefully be chosen for Vab is constrained to a considerable extent by the requirement that this 
equation have solutions for which have the structure exhibited in ( 3.12 ) and ( 3.14 ) for the full 2-point functions 
and inherit the various properties that we discussed in Section III. Observe first that the S{t — t') on the right of (4.2) 
arises from differentiating 9(t — t') and 9(t' — t). We can ensure that only these 5 hmctions will appear by restricting 
V to the form 



V = 



ij°dt 
-i-/°dt 



(4.3) 



where the ellips is indicates terms without time derivatives. The coefficients ±170 are determined by the boundary 
condition ( 3.19 ). In principle, an ansatz using more time derivatives (along with further boundary conditions to 
eliminate unwanted S functions an d th eir derivatives) might be possible, but w e hav e not found such a generalization 
tractable. Next, if the solution of ( |4.2| ) is to have the symmetry expressed by ( 3.15| ), then 2? must satisfy 



which ensures that Iqo is CP invariant. Similarly, the symmetry expressed by ( 3.16| ) implies 



-V22{t,-dt) 

-V2i{t,~dt) 



-Vi2{t,-dt) 
-Vn{t,-dt) 



Finally, as for a scalar field, causality requires 

Pii (i, dt ) + Vu (t, dt) + V21 {t, dt) + V22 {t,dt)^0. 
The net effect of these considerations is that T) can be written as 



V{t, dt 



lV2{t)_ 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



where 'Di{t) and I?2(^) are subject to the constraints 



Pi(t) =I?i(t) 

V2{t)^V2{t) 
Vi{t) + iV2{t) = Vi{t) ~ iV2{t) . 



(4.8) 



When V has this structure, and the propagator S{t,t') is written in the form ( |3.12 ) in terms of a function H{t,t') 
that approximates TL{t,t'), then equations (4.2) for the propagator reduce to 



[ifdt + Vi{t) + iV2it)] H{t, t') = (4.9) 
[i-f°dt + Vi{t)] H{t', t) + tV2{t)H{t', t)=Q. (4.10) 

At this point, the most general procedure would be to expand 'Di{t) in terms of a complete basis of Dirac matrices 
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= ^<(t,fc)rp, (4.11) 

the thirty-two undetermined functions d^{t, k) being the analogues of the functions ak(t), (3k{t) and 7fc(t) that appeared 
in the scalar theory. This general problem is one that we have not found tractable. To simplify matters, we use 
instead the smallest subset of the Dirac algebra that closes under multiplication and under and ~ conjugation, and 
that includes the matrices 7*^ and 7 • k that appear in the free theory. For the operator I?, a convenient basis is 
{l,70,r+,r_}, where 

r± = ^(i±/)7-fe (4.12) 

and 1 denotes the unit matrix. The conjugates of these matrices are 7" — —7'^ ~ 7", V± ~ Tip and T± = T±. Our 
ansatz for T) is then 

V^{t) = I [A(t) - y{t)] 7° + [a{t) - e{t)] T+ + [cj\t) + e*{t)\ r_ - [r(t) + ii^{t)] (4.13) 
V2{t) = ii^{t)j° + e(t)r+ - e*(<)r_ + ir]{t) . (4.14) 

Although we have not indicated it explicitly, the coefficients depend on k and fi as well as on t. For complex functions 
of /i, we define f^{fi) = /*(— ^i), so that f{fi) — f{—fi) = /*''(/i). This is the most general ansatz that satisfies the 
restrictions ( |4.8| ) on Vi (t) , provided that 

a' = a, — T, u* = rj* = r/ 
X + X* =,y + J 

a--(7 = e-? (4-15) 
T — r* = —i{r] — if) . 



V. SOLUTION FOR THE REAL-TIME QUASIPARTICLE PROPAGATORS 



Having constru cted t he unperturbed action (4A) in terms of the differential operator Vj i, dt) given by ( |4.7| ) together 
with the ansatze (4.13) and (4.14), we require a formal solution to equations (^^) and ( 4.10| ) for t he m atrix-valued 
function H{t,t';k) from which the quasiparticle propagator S''°'''\t,t';k) is to be constructed via (3.14). It proves 



convenient to reorganize our basis of Dirac matrices into the set {7± = i(l±7°),r±}, expanding H{t,t') as 
H{t, t') = exp X{t")dt"^ [A{t, t')-f+ + B{t, t')T^ + C{t, t')T+ - D{t, i')7-] 



(5.1) 



Noting that (O) governs the dependence of H{t,t') on its first time argument, while ( 4.10 ) refers to the second time 

o I . 

argument, we introduce the notation A{t,t') = dtA(t,t') and A(t,t') = dt'A{t,t'). With this notation, (4.9) can be 
written as 



(5.2) 
(5.3) 







'A 












'c 









where 



— r 



(5.4) 



For orientation, we note that, in the absence of the counterterm AAab, we would have t = m and a — — |fe|. (Equations 
somewhat analogous to these have been obtained, for example, by Sahni JTot in the course of solving the Dirac equation 
in certain curved spacetimes.) The matrix T has the generalized Hermiticity property — T, where the operation 
is defined by taking the transpose, and replacing each element by its " conjugate. (This reduces to the usual Hermitian 
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conjugate when the chemical potential vanishes.) If the two-component vectors <p and x are solutions to iip = Tip, 
then it is simple to see that the inner product {ip, x) = f'^X is preserved by the time evolution. The eigenvalues of T{t) 
are ±ri(<), where the time-dependent frequency O = Vr^+o^ satisfies $7' = fi, and the corresponding normalized 
and mutually orthogonal eigenvectors are 



1 
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1 



Given some fixed time to, one can formally write exact solutions 



which are positive- and negative-frequency solutions 



T{t')dt' 



to 



-{n-r) 



exp 



Ti / n{t')dt' 



(to), 



at times near to. More generally, we may choose an orthonormal basis 



fHt) 



(5.5) 



(5.6) 



(5.7) 



(5.8) 



with f'^(t)f{t) + g^{t)g[t) = 1. It is readily verified that 4>2{t) is a solution if (f>i{t) is, and that if 01 (i) = ^f,^ "*(/:) for 
some to, then 02(0 = 4'\f, Thus, the solution to (5.2) and (5.3) can be written as 



A{t,t') 
B{t,t') 

C{t,t') 
D{t,t') 



I{t)\ 

9{t) ) 



P2{t') 

'Q2{t') 



-gHt) 
Pit) 



(5.9) 
(5.10) 



where Pi{t') and Qi{t') are to be determined by solving ( I.1C| ) and applying suitable boundary conditions. 
To solve (4.10), it is helpful to define 

W{t,t') = B\t,t') - C{t,t') 
X{t,t') ^ D\t,t') + A{t,t') 
Y{t,t') ^ B\t,t')+C{t,t') 
Z{t,t')^DHt,t')-A{t,t'). 



In terms of these functions, (4.1C ) is 



-i{X + X*) 



T' 



The matrix T is the same as the one defined in ( p.4[ ) , while 



T' 
E 



-e* 



i{v - v"^ - 1] - if) e*f-|-e 

-(e*+ef) i{v - + ri + -q^) 



(5.11) 
(5.12) 



(5.13) 
(5.14) 

(5.15) 
(5.16) 



The appearance of a new matrix T' (which also satisfies T'^ = T') and a new damping constant A + A* in (5.14) 
refiects the fact (previously noted in p2|) that particles and antiparticlcs acquire different thermal masses and widths 
i n the presence of a nonzero chemical potential. The appearance of the inhomogeneous term proportional to E in 
(5.14) implies that the time dependence of the propagators is not exhausted by that of the single-particle mode 
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functions, and we hope to interpret the additional time dependence in terms of time-dependent occupation numbers 

for the single-particle modes^ 

Evidently, the solution of (5.14) will require the introduction of a new set of antiparticle mode functions, which make 
matters rather complicated. In the remainder of this paper, we shall simplify our calculations by specializing to the 
case of zero chemical potential. In th at ca se, " conjugation reduces to complex conjugation, while the restrictions on 
the countcrterm parameters noted in ( ) imply that r, A and rj are real, and that i/ = X. We now have T' = T = T'^ 
and 



-2iri 

-2e* 



2e 
2ir/ 



We look for a solution to ( 5.13 ) and ( 5.14 ) in the form 



(W{t,t') 
\X{t,t') 

(Y{t,t') 
\Z{t.t') 

and find that they are satisfied if 

Si{t,t') 

S2{t,t') 

where 



(f{t') 



\9{t') 
^5.(M')(j[;'j 



+ S2{t, t') ^ f*^t') 



Si{t)L{t') 
S2{t)L{t') 



Ri{t)Ki{t') + R2{t)K2{t') 
Ri{t)K;{t') - R2{t)K^{t') 



L{t')= exp 
m') = Lit') 



\{t")dt" , 
L-\t")e,{t") dt" 



and the ei{t) are related to the matrix elements eij{t) ~ (f)\{t)E{t)(t)j{t) by 



en 
ei2 



2 [e/*.9 - - ■iv{t'){f*f - 9*9)] = iei 

[ej / +^ 9 9 + tVJ 9 \ = . 



(5.17) 

(5.18) 
(5.19) 



(5.20) 
(5.21) 



(5.22) 
(5.23) 



(5.24) 
(5.25) 



At this point, the solution of (O) ( expre ssed by (5^) a nd ( 5.3)) h as lef t us with four undetermined functions 
Pi{t') and Qi{t'), while the solution of ( 4.10| ) (expressed by ( j.lSj ) a nd (5.14 )) produced another four undetermined 
functions Ri{t) and Si{t). However, the two solutions are related by (5.11) and by comparing them, we determine all 
eight functions up to const ants of inte grati on. Moreover, the values that these constants can take are constrained by 
the equal-time conditions ( |3.17 ) and (3.19), which express very general properties of the 2-point functions. In fact, 
it turns out that that the propagator is determined up to two constants of integration, that we shall denote by oi 
(which is real) and a2 (which is complex). Our final result for H{t,t') is expressed by the original ansatz (5.1) with 
the time-dependent coefficients given by 



A{t, t') = [1 - N{t')] f(t)nt') + N{t')g*{t)g{t') - A{t')f{t)g{t') - A* {t')g* {t)nt') 
B{t, t') ^ [1 - N{t')] g{t)nt') - N{t')r{t)g(t') - A{t')g{t)g{t') + A* {t')/* (t)/* (f) 
C{t, t') = [1 - Nit')] f{t)g*{t') - N{t')g*(t)f{t') + A{t')f(t)f{t') - A* {t')g* {t)g* [t') 
D{t, t') = [1 - N{t')] g{t)g*{t') + N{t')r(t)f{t') + A{t')g{t)J{t') + A*{t')r{t)g*{t') 



(5.26) 
(5.27) 
(5.28) 
(5.29) 



The real function N{t') and the complex function A{t') are given, in terms of the quantities introduced above, by 



N{t')^\[l-K,{t')-a,L{f)] 



A{t') = -[K2{f) 



a2L{t') 



(5.30) 
(5.31) 



but it is more illuminating to observe that they obey the differential equations 
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dt'Nk{t')^-2Xk{t')Nk{t') 



(5.32) 
(5.33) 



with initial conditions A^/c(0) = ^(1 — aik) and Afc(O) = ^a2k, in which we have reinstated the dependence on spatial 
momentum k that has been implicit throughout. These are the equations that we might hope to interpret as kinetic 
equations for the occupation numbers of single-particle modes. Let us, indeed, specialize to the case of thermal 
equilibrium, and temporarily delete the counterterm M.ab- The mode functions can be written as 



fk{t) 
9k{t) 



(5.34) 



with Qk = Vk'^ + ra^ and we have itjfe = K2k = and Lk = 1. We then find that that the time-ordered function 
S'("'(t,i', fe) = H{t,t'J{)9{t - t') + H{t',t;k)e{t' - t), with H{t,t';k) given by (|1|) agrees with the corresponding 
function obtained in |15t] , provided that we can identify Nk = ^(1 — ai) 



and Ai. 



J 02 = 0, where N. 



[exp(/3f2fe) + 1] is the usual Fermi-Dirac distribution. The other 2-point functions do not agree with those of [|5|, 
because these authors made use of a different time path (which is legitimate in thermal equilibrium, but not in the 
nonequilibrium situation considered here). At this stage, these values of the constants of integration oi and 02 are 
merely guesses that yield this agreement with . The actual values that are required by our formalism are determined 
by computing imaginary-time correlators and applying appropriate boundary conditions. This computation is the 
subject of the following section, where we shall find our guesses confirmed. 



VI. IMAGINARY AND MIXED-TIME PROPAGATORS 



In terms of the imaginary-time field operator tp{x,T) — e^'^il){x,0)e , the imaginary- and mixed-time 2-point 
functions are (for a — 1,2) 



' (x, t; x', t') = {i,^{x, T)i;p{x\ r'))0(r - r') - {^^{x', r')Va(a;, T))0(r' - r) 
S^^^\x,t-x\T') = -{^p{x\r')Mx,t)) 
5(3") (a;, r; x', t') = S^"^'^ (-a;', t'- -x, t) . 



(6.1) 
(6.2) 
(6.3) 



Because the time-path ordering makes imaginary times later than real times, the functions S'^^'^^ and iS^^"^) are identical. 
With /i = 0, antiperiodicity of the path integration variables, ?/'3(a;, p) — —ipi{x, 0), and the fact that the field operator 
'ijj{x^ 0) is unique supply the two boundary conditions 



5("3)(t,0;fc) =5("2)(t,0;fc) 



(6.4) 
(6.5) 



that we shall use to determine the constants of integration ai and 02. If the sources for real-time fields are set to 
zero, then the time path reduces to just its imaginary-time segment, and antiperiodicity yields 



5(33)(0,r';fc) = -5(33)(/?,r';fc). 
Finally, uniqueness of ipi^^ 0) ^-Iso implies 

5('^3)(o,r';fc) =5(33)(0,r';fe) 



(6.6) 



(6.7) 



and the two latter boundary conditions serve to fix constants of integration that arise in the calculation of 5*^"^) and 

In order to construct a tractable perturbation theory, we have insisted that the unperturbed action (4.1) be local 
in time. With this restriction, there are no terms of the form ^ ^^(tyD a^ipziT) , {a = 1, 2) so we have 



/co(V^) 



E 

a, 6=1 



dt tp^Vabipb + dr V's^'saV's 



(6.8) 
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In the second term, which approximates the path integral representation of the initial density operator, ip^'D^^ip^ is 
—i times the Euclidean version of the free part of ( ^.l| ), with m = m(0), supplemented by a counterterm A^sa. The 
form of I?33 is determined by the assumed CP invariance, I?33(9t) — 2^33 (—9-r), together with the requirement that 
the quasiparticle energy Qq be equal to the t limit of the energy that appears in the real-time mode functions. 
This yields 



Vssidr) = I [l^dr - (T^T+ - C7o*r_ + Tq] , 



where ctq and tq are the t ^ Q limits of the parameters appearing in (4.13). The new propagators satisfy 



(6.9) 

(6.10) 
(6.11) 



and these equations can be solved by the method explained in the previous section. The imaginary-time propagator 
can be expressed as S^"^^^ (r, r'; fc) = H^{t, t'; k)9{T — r') + H^{t' , r; k)6{T' — r), and we look for solutions of the form 



i/3(T, r') = Asir, t')7+ + Bsir, r')r_ + C^ir, t')T+ - D^ir, r')7- 
^(13) (i^/) ^ exp (^-^ \it")dt"^ [Ai3(i,T')7+ +Si3(t,r')r- +Ci3(i,r')r+ - Dr3{t,r')j_] 

Defining positive- and negative-frequency mode functions in imaginary time by 



Mr) 
9i{r) 



-Hot 



Mr) 



V2rio(r!o - To) 



'(f^o - To) 



the solutions subject to the boundary conditions ( |6.6D and (6.7) are given by 

A3(r, r') = (1 - iV^'^)/i(r)/i(r') + N'^^g,{T)gi{T') 
B3{r,r') = (1 - N^^)g,{r)Mr') - A^'=Vi(r)gi(r') 
C3(r,r') = (1 - 7V^'^)/i(T)gi(r') - 7V'='igi(r)/i(T') 
i?3(r, r') = (1 - N^^)g,{T)g,{r') + iV^Vi(r)/i(T') 



and 



yli3(i,r') = -Ar°q/(t)/j(r') - (1 - iV^'i)5*(t)<7i(r') 
i3i3(t,T') = -N-^g{t)Mr') + (1 - A^°^)/*(0.9i(t') 
Ci3(i,r') - -7V°V(i)5i(r') + (1 - A^'='').9*(t)/i(r') 
i^i3(i,T') = -7V°^g(i)5i(T') - (1 - iV'^'^)/* (0/i(t') , 



(6.12) 
(6.13) 

(6.14) 



(6.15) 
(6.16) 
(6.17) 
(6.18) 



(6.19) 
(6.20) 
(6.21) 
(6.22) 



with N'^'i = [e^"« 
ai = 1 - 2N°'i and 02 



[e^^" + 1] ^. Finally, the boundary conditions (6^) and (|6.5| ) are satisfied provided, as promised, that 



VII. CONCRETE REALIZATION: A SIMPLE MODEL 



In the preceding sections, we have constructed a quasiparticle action, and the propagators that correspond to it, in 
terms of several time- and momentum-dependent coefficients that so far are undetermined. In general terms, these 
coefficients are to be determined self-consistently by asking the counterterm Aiab to cancel some part of the higher- 
order contributions to the self energy, which arise within an interacting theory that we also left unspecified. Thus, 
the full 2-point functions can be expressed through the Schwinger-Dyson equation 

S'^''''\t,t';k) = S'^''''Ht,t',k)~i J dt"dt'"S^'"'Ht,t";k)j:ad{t",t'";k)S^'^''\t"',t';k), (7.1) 

in terms of a self energy that has contributions both from the counterterm and from loop corrections 

Sab(i, t') = Mabit, dt)6{t - t') + ^'°°P{t, t') . (7.2) 
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If the counterterm could be chosen so that Sab = 0, then the propagator S'('^f') would be the same as the full 2-point 
function 5^"''-' and the perturbation series would be completely resummed. In practice, of course, we can achieve 
at best a partial resummation by cancelling the first few terms in the expansion of S'°°p. Moreover, since T,^°°p is 
non-local in time, and may well have a more complicated spinor structure than the counterterm we have constructed, 
it will not be possible to cancel even these terms exactly. We can effect a selective resummation by cancelling only part 
of E'°°p, but the choice of which part to cancel will depend on details of a specific application and of the supplementary 
approximations that will inevitably be required. Here, we illustrate how the process can be made to yield sensible 
results by studying the simplest possible model, in which our fermion interacts with a real scalar field, the interaction 
being specified by 



(7.3) 



We suppose that the (f> particles have a mass M that is greater than 2m, so that the decay and annihilation processes 
(j) ^ i/; + i/j are kinematically allowed. We anticipate that these on-shell proces ses w ill give absorptive parts to the 
fermion self energies, yielding a nonzero thermal width A, and that the equation ( 5.32| ), which gives the rate of change 
of the occupation numbers Nk{t), will, within a suitable approximation, be recognisable as a kinetic equation of the 
Boltzmann type. 

We study explicitly the one-loop, real-time self energy, corresponding to the emission and reabsorption of a 
particle. By setting M + I]1-i°°p « 0, we obtain a complicated set of constraints, which implicitly specify the 
functions Xk{t), ek{t), etc. In fact, these functions enter 5]^~'°°p through the mode functions fk{t) and gk{t) and the 
auxiliary functions Nkjt ) and Afc(^ ) , for which we have no concrete expressions in hand. We know only that they are 
solutions of (pT^), (5.35) and (5.3S). In principle, we have a closed set of equations that we might attempt to solve 
numerically. To see more clearly what these equations imply, however, we introduce some further approximations. 
First, we will suppose that time evolution is sufficiently slow for an adiabatic approximation to be reasonable. Then 
the mode functions will be written as 



fk{t) 



1 



9k(t)J ~ ^2nk{t){nk{t) ~ Tk{t)) \nk{t)-Tkit) 



<Jk{t) 



exp 



i / Vlk{t')dt' 



(7.4) 



Further, when the fermions have a nonzero thermal width, the propagator will decay at large time separations, very 
roughly as e""^'*"* L Assuming that Tk{t) and ak{t) do not change to o mu ch over a thermal lifetime \k{t)~^, it will 
be reasonable to approximate the product f{t)f*{t') that appears in ( 5.26 ) as 



ak{t)<Jl{t) exp [-ink{t){t - t')] 
2^k{t) [^k{t ) - Tk{i)] 



(7.5) 



where t — {t-\- 1')/2, with corresponding approximations for other products of mode functions. Finally, since X1^^'°°p 
is proportional to , so are the functions Afc(t), etc that appear in A^. At the lowest order of perturbation theory, it 
is therefore reasonably consistent to set these functions to zero in the propagators that we use in evaluating 5]i-ioop^ 
and this is what we do. In particular, we then have ak(t) = —\k\ and Tk{t) = m{t). Clearly, these approximations are 
valid, at best, only for a weakly interacting system in a state close to thermal equilibrium. It is therefore important to 
emphasise that our purpose in introducing them is to obtain simple analytical results that illustrate essential features 
of our formalism. The formalism itself is by no means restricted to situations where these approximations are valid. 
For different reasons, discussed below, we will set A(i') — 0. 



At this level of approximation, the propagators used in evaluating 1]1-i°°p (and those that multiply Safe in (7.1)) are 
essentially those of the equilibrium theory, except that we allow for time-dependent masses and occupation numbers. 
It is useful to introduce the projection operators 



A(fc) = r2fe(t )7° - 7 • fc + m{t ) 
A(fc) = -rjfc(t )7° - 7 • fc + m{i) 
'm'^it) - P 



A±(fc) =± 



7" - 7 • fe + m{t) 



(7.6) 
(7.7) 

(7.8) 



which have the properties 



A(fc)A(-fc)A(fc) 
A(fc)A(-fc)A(fc) 



A(fc)A_(A;)A(fe) = 
A(fc)A+(fc)A(fe) = 0. 



(7.9) 
(7.10) 



12 



After a Fourier transform on the time difference t — t', each of the propagators appearing in the second term of (7.1) 
can be written in the form 



(fe, w; ? ) « si"'"'' {k,uj; t)A{k) + s!f''^ {k,uj; i ) A(fe) 

(ab). 



(7.11) 

where s[°'^\k,uj;i) has poles at w = flk{t)±iXk{t), while S2"'''\k,u!;t) has poles at u; = —nk{t)±iXk{t) and, within 
the approximations described above, Xk{t) is to be regarded as infinitesimal. In particular, we shall make explicit use 
of S^^^^k, uj; t ), in which the poles combine to yield 

(7.12) 



(7.13) 



^ 1 [^^0 + m{t)] [Nkit)S (uj ~ 17fc(i )) + (l - iVfe(f)] s{lo + ilkit)^ 



and of the corresponding scalar propagator 



UJk{t) 



nk 



{t )S (lu - uok{t )) + (l + "fc(f )) 5{uj + ujk{t) 



with Wfc(0 = y/W+MHt). 

We now consider how our vaguely-stated criterion + I]'°°p k, can be implemented in practice, to yield a selective 
resummation of the perturbation series. To be clear, let us first reiterate the role of the simplifying approximations 
introduced above. In principle, the functions Tk{t), <Tk{t), Xk{t) ■ ■ ■ are to be determined self-consistently by solving a 
set of equations of the form A4(t, ct, A, . . .) w — I]'°°p(t, a,X, . . .). However, for the purposes of discovering how « might 
sensibly be interpreted and eventually of making contact with kinetic theory, we plan instead to study the simplified set 
of equations A1(t, ct, A, . . .) ~ — I]'°°P(m, — |fc|, 0, . . .), for which we can find closed-form expressions. Of the functions 
that we wish t o d etermine, r and a clearly encode the dispersion relation for the quasiparticl es wh ose mode functions 
are given by ( [5.6[ ). These functions appear only in Vi (for which our ansatz was given in (4.13)) and thus in Aiu 
and 7Vi22- By arranging for them to cancel appropriate parts of 1]^^°^ and S2°2°''j endow our quasiparticles 

with dispersion relations that approximate those of the true elementary excitations of the nonequilibrium state. We 
do not do this explicitly, however, preferring to focus on dissipative aspects of the problem, and especially on the 
evolution of occupation numbers. In fact, for simplicity, we shall set Tk jt) =_m{t) and ak{t) = — |A;| in M. as well 
as in S'°°P. The terms proportional to Afc(t) in our propagators ( f3.26| ) - ( |5.29| ) are awkward and (as discussed in 
the next section) difficult to interpret. Within our present approximations, it is possible to eliminate these terms in 
the following way. We saw at the end of section VI that the constant a2k = Afc(O) vanishes when the initial state 
is one of thermal equilibrium. It is therefore consi stent to set Afc(t) = at all times, provided that the quantity 
&2k{t) vanishes in (5.33). The expression given in (5.25) does not vanish in general, but with the approximations 
made here it can be made to do so by choosing ek{t) to be purely imaginary, say ek{t) = iik{t), and by choosing 
Vk{t) = -Tk{t)ek{t)/<Jk{t) = m{t)ek{i)/\k\. 

With these simplifications, the functions that remain to be determined are Xkit) and ikit), which appear in the 
counterterm 



Xl2(fc,0 



2^k{t) 



Afe(t) + ^^.W) A(-fe) - (Xk{t) - ^e,(t) ) A(~fe) 



(7.14) 



Ideally, we would like this to cancel the one-loop contribution 



t1— loop 
-•12 



{k,uj;t) w -ig^ 



duj'duj" fd^k'd^k" 



27r 



(27r)3 



S{uj -uj'- u")5{k k' k") 5(i2)(fc', uj'- t)g'^^^^ {k", uj"- i) , (7.15) 



but M.12, being derived from a counterterm that is local in time, has no dependence on w, while cannot be 



expressed in the form of ( 7.14 ) as a linear combination of A(— fc) and A(— fe). Clearly, we must be a little less ambitious. 
First, we shall attempt to effect the desired cancellation on shell: that is, to cancel only sjg^'""'^ {k^±D.kit)^ . For 



Lo = ±r2, we find 

^1— loop / 



E™(fc,c.;f) 



167r2 



+ ^7^7 ■k + m]n{l- N')9{uj) 



k ■ k' 



(7.16) 



Here, the kinematics is that of on-shell decay or pair annihilation, involving two fermions of momenta k and fc', with 
energies J7 = ^lk{t) and H.' — rife'(f), and a scalar with momentum p = k + k' and energy ujp = + M^{t). The 
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fermion occupation numbers N = Nk{t), and A'"' = Nk'{t) are those that we originally introduced in (5.26) - (5.29) 
while n — Upit) is the corresponding quantity for the scalar 



When uj is close to +f2fe, the propagators given approximately by (7.11) can be further reduced by retaining only 
the term containing a pole, namely S'(°^) « S'j"''-*A(fc). The one-loop correction term in the Schwinger-Dyson equation 



(7.1) then involves the projection A(fc) A^i2 + A(fe). Using the properties (|7.9| ), we see that this can be 



made to vanish by expressing eJ^'°°'' as a linear combination of A(fc), A(— fc) and A_(fe) and requiring the coefficients 
of A(fc) in ^'^d M.12 to cancel. Similarly, when lo is close to —fife, we express eJ^'°°'' as a linear combination 

of A(— fc), A(— fc) and A-|_(fe), and require the coefficients of A(— fc) to cancel. In this way, we obtain 



647r2 ' ' J nn'LUp 

^2 



(7J7) 



Here, we have made use of the kinematic identity QQ' — k ■ k' ~ w? — ^(M^ — Aw?). Reassuringly, we have arrived 
at a thermal quasiparticle width that is positive-definite when M > 2m, so that the on-shell decay and annihilation 



processes that give rise to it are kinematically allowed. When M < 2m, the thermal width v anish es, because (7.15) 



contains products of (5-functions that cannot be satisfied simultaneously. The expression (7.17) agrees with the 
equilibrium damping rate calculated in ||l^ (see also p8j)for t he sa me model. 



Finally, we can evaluate the right hand side of equation ( ^.32 ), which we hoped to interpret as governing the 



evolution of time-dependent occupation numbers Nj^it). With the approximations used in this section, we have 
eife(i) — — 2f2fe(t)efc(t)/|fc|, and the equation becomes 

dNkit) _ g2 ^^^^^ _ ^^^^ r ^^^^ S{ujp~n-n') ^^^^ _ ^^^^ _ _ ^ ^^^^^^ 



dt 327r2' 'J nn'ujp 

We recognise the standard form of a relativistic Boltzmann equation (restricted to a spatially homogeneous system) 
in which the gain and loss terms have the correct statistical factors, n{l — N){1 — N') and (1 -|- n)NN' respectively, 
to represent a fermion of momentum fc being produced by the decay of a in the thermal bath or annihilating with 
an antifermion in the bath. Again, the 'collision' integral would be replaced by if M < 2m. Of course, our results 
for Afc(t) and dNkit) /dt were obtained only at the lowest non-trivial order of perturbation theory. At higher orders, 
we would expect non-zero answers in both cases, arising from scattering processes that are allowed even for M < 2m. 

VIII. DISCUSSION 

We have described a selective resummation of the perturbation series for the nonequilibrium 2-point functions 
of spin-i fermions. The general philosophy of this resummation is to describe the nonequilibrium state as nearly 
as possible in terms of its own quasiparticle excitations. The propagators for these excitations, unlike the free- 
particle propagators used in standard perturbation theory ought, roughly speaking, to incorporate nonzero widths 
and occupation numbers that evolve with time, reflecting the evolution of the nonequilibrium state. The resummation 
is achieved through the use of a counterterm that transfers some contributions of higher-order self energies into the 
lowest-order theory about which we perturb. In a nonequilibrium situation, this is tractable only if the unperturbed 
action is local in time, and this places strong constraints on what can be resummed in practice. For example, 
resummations somewhat similar in spirit to ours, but restricted to scalar theories in thermal equilibrium, are described 
in [^-|T] and applied to the 'warm inflation' scenario in In thermal equilibrium, the 2-point functions depend 
only ont — t', and after a Fourier transformation, one can (in principle) construct a counterterm analogous to M that 
subtracts the whole frequency-dependent self energy at whatever order of pe rtur bation theo ry one has the energy to 
compute. Generalizing this to a nonequilibrium state would mean replacing (U) and ( |4T0| ) with integro-differential 



equations containing arbitrary non-local kernels in place of the functions Tk{t), Xkit), etc. (or perhaps with infinite- 
order differential equations having infinitely many time-dependent coefficients, all to be determined self-consistently) 
and this does not seem to be a practical proposition. For fermions, this may be particularly unfortunate, because 
the high-temperature plasma has, at least in some important cases, 'hole' or 'plasmino' excitations in addition to 
the particle and antiparticle poles that are present at T 0. These have been known for some time from one-loop 
calculations ||^ and their properties have recently been explored in terms independent of perturbation theory by 
Weldon p^ ] . The counterterm we have constructed is linear in dt and cannot accommodate a dispersion relation with 
these multiple branches (though a generalization that mimics them might be possible). 



14 



Despite this inevitable deficiency, the resum mati on appears to make reasonable sense. Of the initially arbitrary 
parameters that we introduced in ( 4.13| ) and ( 4.14| ) to describe the quasiparticle excitations, Tk{t) and ak{t) have 
clear interpretations in terms of the quasiparticle dispersion relation, while Xk{t) is a thermal width which, after 
enough weak-coupling and adiabatic approximations, turns out to agree with the fermion damping rate calculated 
in equilibrium. In the same approximation, the function Nk{t) that appears in our resummcd propagators does 
indeed correspond to a quasiparticle occupation number, evolving according to a kinetic equation of the Boltzmann 
type. It is worth emphasizing that, while this kinetic equation provides reassurance that our formalism has a sensible 
interpretation, its derivation is by no means the purpose of the formalism presented here. There are indeed, many 
routes to equations of this kind. One may, for example, investigate directly the time evolution of the expectation value 
of a time-dependent number operator (see, e.g. ||2^ ). Another route that has been pursued extensively in connection 
with the calculation of transport coefficients of high-temperature plasmas is to extract a transport equation for a 
Wigner density through truncation and gradient expansion of the Schwinger-Dyson equations |25-28|. Attempts to 
calculate transport coefficients directly from the Kubo formula of linear response theory reveal, on the other hand, 
infrared singularities that require the resummation of larg e classes of diagrams [ p9|j30| ] and this resummation turns 
out to be equivalent to solving a Boltzmann equation [pl[ . (The relationship between these approaches is discussed 
in [^2|). The transport equations that arise in these calculations go well beyond the simple one- loop approximation 
exhibited in ( 7.1£ ), but they apply to systems very close to equilibrium. Our own goal of constructing a resummed 
perturbation theory to describe the evolution of highly-excited states that may be far from equilibrium is rather 
different. The functions Nf^{t) that arise in the course of solving for the resummed propagators arc not necessarily 
equivalent to the Wigner distribution, and the equation ( ^.32 ) that defines them reduces to a Boltzmann equation 
only after approximations that one may in general hope to avoid. 

Finally, the investigation reported here indicates that the structure of nonequilibrium fermion propagators is more 
complicated than might be expected from the equilibrium theory. Close to equilibrium, we found that their spinor 
structure can be expressed in terms of the two projection operators (7.6) and (7.7), as can be deduced on general 

I])- 



grounds within the equilibrium theory (see, e.g. |2^]). Away from equilibrium, this is no longer true. In general, 
there are at least terms proportional to ^'^f ■ k and possibly other terms that we have not succeeded in resumming. 
To arrive at propagators that contain only A(fc) and A(fc), we need the coefficients B{t,t') and C(t,t') in ( pT| ) to 
be equal. This will be true of the solutions presented in (5.27) and (5.28) if the function A{t) vanishes, and if the 
products of mode functions g{t)f*{t') and f{t)g*{t') are equal. Both of these conditions would be automatic if we 
were to assume time-translation invari ance , so that the propagators depend only on t — t' , w hich is, of course, true 
in thermal equlibrium. The equation (5.33) satisfied by A(t) is superficially similar to (5.32), but does not bear the 
same interpretation as a kinetic equation. In fact, the terms in the propagator that involve this function are products 
(loosely speaking) of two positive-frequency or two negative-frequency mode functions that have no counterpart in 
equilibrium. In general, it seems that such terms should be present in nonequilibrium self energies, but we have found 
no simple interpretation for them. 
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